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Analyses of a modified BCS (MBCS) theory performance at finite tem- 
peratures in the Picket Fence Model (PFM) for light and heavy systems 
Ph ! is presented. Both symmetric, = N (N particles on f2 twice-degenerate 

levels), and asymmetric, Q ^ N, versions of the PFM are considered. 
The quantities known exactly from particle-hole symmetry of the f2 = N 
PFM are calculated. Starting from very low temperatures, these quan- 
tities are found in dramatic deviation from the exact values in MBCS 
results. Consequences of the MBCS prediction that heating generates a 
! thermal constituent of the pairing gap, are discussed. Thermodynamical 

inconsistency of the MBCS is also addressed. 
PACS: 21.60.-n, 24.10.Pa, 24.60.-k, 24.60.Ky 
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O ' 1 Introduction 

o , 

A modified BCS (MBCS) theory for treating pairing correlations in atomic nuclei at finite 
temperatures [T|2] has been recently tested [3] in the Picket Fence Model (PFM) in which 
N particles are distributed over Q twice-degenerate levels. The PFM with N = Q is usually 
considered in the literature. 

The MBCS predicts a smooth decreasing behavior for a pairing gap as temperature T increases 
up to some T M when the theory suddenly breaks down. It was reported [3|3] that adding 
one extra level (O = N + 1) extends the MBCS applicability to much higher temperatures. 
As is already established [5], the Q = N + 1 case of the PFM is the only example of an 
exceptional MBCS pairing gap behavior and thus, it could be hardly considered as a typical 
one. Accordingly, we find it necessary to provide the reader with other Q ^ N PFM examples 
not available in [3]. This is done in Sec. [21 

In Sec. El the MBCS prediction that not only the pairing force but also the heating itself 
generates the pairing gap, is discussed. 

In Sec. HI we address the key question of the article: whether the MBCS is a reliable theory 
in the temperature domain where its pairing gap looks reasonable at first glance. For that, we 
examine some quantities which are known exactly in the N = Q PFM because of symmetry. 
We show that in the MBCS predictions, these quantities dramatically deviate from their exact 
values starting from very low temperatures. 

Thermodynamical inconsistency of the MBCS is discussed in Sec. [51 



2 MBCS pairing gap in PFM systems with Q = N and Q ^ N 



The PFM or Richardson model is widely used as a test model for the pairing problem. It is the 
pairing Hamiltonian applied to a system of N fermions distributed over Q equidistant levels. All 
levels are twice degenerate for the spin up and down. The levels below (above) the Fermi surface 
will be referred to as holes (particles) and labeled by "— i" (V). Their single particle energies 
are e_j = (0.5 — i) MeV for holes and £j = (—0.5 + i) MeV for particles, where i — 1, 2, , . . . 
(i.e. Si = —£-{)■ In all calculations presented below (except for the ones in Fig. E^), the pairing 
strength parameter G is adjusted so that the pairing gap A equals 1 MeV at zero temperature. 

In addition to the fact that exact solution of the PFM is possible if N and Q are not big, the 
PFM with N = Q (to be referred to as the conventional PFM) possesses internal particle-hole 
symmetry. This means that at any temperature: 

a) the energy of the Fermi surface Ep equals exactly MeV: 



Ep = 



(1) 



b) the quasiparticle energies Ei for particles and holes should be degenerate: 

Ei = 
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because of (CQ); 

c) the particle occupation probabilities (Bogoliubov Ui and Vi coefficients) are related as: 



Ui 



E, 



E, 



\ 




V-i 



(3) 



because of ( TTf2"|) ; 

d) the thermal quasiparticle occupation numbers = 1/(1 + exp (Ei/T)) and quasiparticle- 
number fluctuations SNi = \Jriiil — rii) should be equal for particles and holes with the same 
i: 

n, = n_i and 5Ni = 5N_i (4) 



because of (|2j). 

Of course, an asymmetric version of the PFM with Q = N + k (where k = —N/2 + 1, . . . , — 1, 
1, . . . , oo) may be considered as well but Eqs. (PH]) are not valid for it. 

The first test of the MBCS performance in the conventional PFM with N = Q = 10 revealed 
that at T ~ 1.75 MeV the system undertakes a phase transition which manifests itself as a 
sharp simultaneous increase in the pairing gap, a sharp decrease of the system energy and a 
discontinuity in the specific heat Cy (this phase transition was defined as a superfluid - super- 
superfluid phase transition) [6] . This critical temperature was denoted as T M [5]T4] . It was found 
that Tm linearly increases with the number of particles N in the conventional PFM. 

It was also reported that enlarging the space by one more level, £1 = N + 1, restores the MBCS 
applicability to much higher temperatures even for iV < 14 systems [5|4] . Fig.[T]plots the MBCS 
pairing gap for iV = 14 particles and Q = N + k levels with k changing from —4 to 50 (solid 
curves). Indeed, the MBCS gap above T c is small and almost constant up to Tm = 7.2 MeV 




12 3 

T (MeV) 



Fig. 1. MBCS pairing gap A for the system of N = 14 particles on 0, = N + k levels (solid curves). 
The BCS pairing gap is plotted by the dashed curve. The k values are indicated in an oval at each 
curve. The conventional critical temperature T c is shown in this and all other figures by the vertical 
arrow. 

when suggested example of extended configuration spaces, k — 1, is considered. However, the 
MBCS theory breaks down again at rather modest temperature Tm = 2.15 MeV with adding 
one more extra level, k = 2, etc. Thus, the extension of the configuration space beyond k = 1 
makes the MBCS results inappropriate again. 

One immediately notices from Fig. [T] that the MBCS theory predicts two typical scenarios for 
the system evolution with heating. As temperature increases, the system undertakes either 

a) a superfluid - super-superfluid phase transition (examples with k < in Fig. [[]) 
or 

b) a phase transition from a superfluid phase with a positive gap to another superfluid phase 
but with a negative gap (examples with k > 1 in Fig. CD), 

instead of a superfluid - normal phase transition of the conventional BCS. Corresponding BCS 
pairing gap behavior in the same systems is shown in Fig. [T]by the dashed curve. It is impossible 
to visually distinguish the results when k value changes from —4 to 50. 

We have performed additional calculations for the asymmetric PFM with N changing from 6 
to 100. The results look very similar to the ones in Fig. [TJ In all these examples, there exists 
only a single case k — 1 with abnormally large Tm which grows almost linearly with N. For 
all other k values, the phase transitions of unknown types take place at much lower Tm which 
becomes smaller and smaller as k differs from 1 larger and larger. 

However, only the case Q = N + 1 has been selected for presentation in [3114"] as an example of 
an extended configuration space (Q > N). 
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(a)N=Q. = 10 



(c) N=£l = 30 





(b) Af = 10, Q= 11 



(d) N=Cl= 100 
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Fig. 2. MBCS predictions for the ratios of quasiparticle-number fluctuations 5Ni/5N-i in different 
(iV, Q) systems. Exact result 5Ni/5N_i = 1 in Figs. [2]Ja,c,d) is shown by the dotted line. 

3 MBCS pairing gap: quantal and thermal constituents 

Although the MBCS performance has been already discussed in several articles, a clear answer 
to the question what is responsible for phase transitions of unknown types predicted by this 
theory is not yet established. According to the MBCS founders, there exists a criterion of the 
MBCS applicability, according to which the line shape of the quasiparticle-number fluctuations 
5Ni should be symmetric [3]H] . And the theory reaches the limit of its applicability at Tm- The 
physical origin of this criterion is unclear. The physical spectra are never symmetric, at least 
in nuclear physics. On the other hand, the single particle spectrum of the conventional PFM 
is ideally symmetric with respect to a chemical potential, see Eqs. (121141 . but the MBCS itself 
breaks down this symmetry (see below). 

The reader of [3] is suggested to notice from Fig. 2 that the line shape becomes rather asymmet- 
ric approaching Tm, thus violating the above-mentioned criterion of the MBCS applicability. 
We find it difficult to judge of the symmetry of lines in Figs. 2(a-e). On the other hand, it is 
very easy to quantify the results in this figure by plotting the 5Ni/5N_i ratios as a function of 
temperature for different i. This is done in our Fig. [2] for different systems. 

Indeed, the asymmetry increases above Tm in the N = Q = 10 system, Fig. [2]^a), while this 
does not in more heavy symmetric N = Q, Figs. [2](c,d), or selected asymmetric Q — N + 1, 
Figs. [21(b), systems. However, it is difficult to miss that this effect is marginal compared to 
asymmetry in all systems (symmetric, asymmetric, light, and heavy) at much lower T with the 
maximum around T c . In other words, the most severe violation of the suggested criterion of the 
MBCS applicability takes place around the critical temperature of the conventional BCS and 
not at T M where the MBCS breaks down, as stated in [3"f4"] . This conclusion does not depend 
on details of the PFM system being considered. 

Since the suggested criterion does not help to understand what causes the theory breaking 
down, we continue our analysis. Let us read once again the paragraph containing Eq. (26) in 
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Fig. 3. MBCS pairing gap (thick solid curve) for the PFM systems with (a) N = 14, 0, = 15 and (b) 
N = Q = 100 and in calculation with a realistic single particle spectrum for 120 Sn: (c) neutron and (d) 
proton systems. This gap is the sum of a quantal (particles plus holes - thin solid curve) and thermal 
(holes - dashed curve, particles - dot-dashed curve) parts. 

Ref. [3J that an important feature of the MBCS theory is that the MBCS gap is the sum of a 
quantal part (which looks the same as in the conventional BCS) and a thermal part SA: 



holes 

5A h = GY,( v i -u 2 i)5Ni for holes (5) 

i 

particles 

SA P = G Yl ( v i ~ u 1) 5N i for particles . (6) 



Fig. [3] presents the MBCS gap for (a) light and (b) heavy PFM systems by the thick solid 
curve. The quantal part of it is shown by the thin solid curve, it becomes very small soon above 
T c . The thermal part for holes (particles) is plotted by the dashed (dot-dashed) curve. The 
quantity 8 Ah is always positive (vi > v*) and the quantity 5A P is always negative (ui > vi); 
their absolute values increase with temperature. 

The physics of pairing suggested by the conventional BCS is very simple and transparent. 
Pairing is generated by the pairing force. As temperature increases the thermal scattering of 
nucleons becomes stronger and stronger and finally destroys the pairing at T c . 

The MBCS suggests another physics, according to which the heating itself generates a thermal 
constituent of the pairing gap: positive for holes and negative for particles. The heavier is the 
system, the stronger thermal pairing gap may be generated. A similar phenomenon takes place 
in calculation with realistic single particle spectra (see Figs. [3(c) and (d) where the pairing gap 
behavior in 120 Sn is presented for neutrons and protons, respectively). 

Notice, when the pairing strength is weak to generate pairing at zero temperature, the MBCS 
predicts that the heating develops the pairing gap at finite T, as is obvious from Eqs. ([MS]). An 
example of the pairing induced by heating alone in a magic nuclear system is shown in Fig. [3jd). 
One should not be surprised that the resulting pairing gap is negative in this example. Fig. [3] 



clearly shows that as temperature increases, the MBCS gap receives main contribution from 
5 Ah and SA P terms and their sum can be equally positive (less particle levels, k < in Fig. [I]) 
or negative (more particle levels, k > 1 in Fig. [[]). 

The main goal of the MBCS theory is to mimic the thermal behavior of the pairing gap of 
a macroscopic treatment [H] that the normal - superffuid phase transition is washed out and 
soon above T c the gap remains rather small but positively finite. However, Fig. [3] shows that 
even technically, it is very difficult to achieve the desirable goal with the MBCS equations: 
it is necessary that two almost linearly growing functions 5Ah and — 5A P almost cancel each 
other with a high accuracy in a large temperature interval^} Of course, the final result is very 
sensitive to tiny details of a single particle spectrum employed. If a spectrum with desirable 
properties is occasionally found, the theory breaks down anyway with a step aside from it. 
Fig. [TJ clearly demonstrates it. 



4 On range of the MBCS validity 

As is established, the MBCS theory is capable of generating the pairing gap behavior which 
looks reasonable at first glance below Tm- It is also established that Tm grows linearly with 
the mass of the system and it is possible to find a single exceptional case of abnormally large 
T M . Does this mean that the MBCS theory should be considered as reliable for T < T M and 
recommended for applications? 

To answer this question we return to the conventional PFM. Due to its internal particle - 
hole symmetry, the accuracy of the MBCS predictions for the theory variables can be easily 
examined even without having exact solutions. 

Verification of Eq. (TJJ for the SNi quantities in Figs. [2](a,c,d) reveals that this property of the 
system is enormously violated in the MBCS predictions. The behavior of the ratios nj/n_j is 
very similar to the one of 5Ni/5N_i in Fig. El Deviation from the exact result nj/n_j = 1 
reaches a few thousand per cent near T c . 

The accuracy of the MBCS predictions for the quasiparticle energies, Eq. ([2]), and for the particle 
occupation probability, Eq. ([3]), with i = 1 is examined in Fig. 0] and Fig. G^a), respectively. 
Fig. E^b) shows that the MBCS theory does not keep the energy of the Fermi surface at zero 
energy, as it should be, Eq. ([1]). 

The results in Figs. [2]^a,c,d), HJ and [5] obviously demonstrate that as soon as heating starts to 
play a role, the MBCS predictions are irrelevant to a system which the theory aims to describe, 
because genuine properties of the system are severely violated in these predictions. It does not 
matter whether a light or heavy systems is considered, the results are almost independent on 
N. 

All-in-all, and this is the main conclusion of the article: although the MBCS theory generates 
a pairing gap which looks reasonable in some temperature interval [0, Tm] , it is clear from the 
analysis of other quantities that the theory predictions have nothing to do with the properties 
of a system under discussion even in this temperature interval. 



1 Independently from the physical content of the thermal gap in the MBCS, one finds from Eqs. ([3ED 
that <5A/j = — 5A p for the conventional PFM and the cancellation between these two terms should 
be exact at any T. This does not happen in MBCS calculations because the MBCS theory violates 
Eqs. (|3|4|) . see below. 




Fig. 4. MBCS predictions for the ratios of quasiparticle energies [E-ijEi) in different N = systems. 
Exact results [E-ijEi) = 1 are shown by the dotted line. 
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Fig. 5. MBCS predictions for (a) ratios [v-\ju\) and (b) the energy of the Fermi level Ep in N = £1 = 
10, 30, and 100 systems. Exact results (a) (v-i/ux) = 1 and (b) Ep = for these systems are shown 
by the dotted line. 



5 Thermodynamic inconsistency of the MBCS theory 



We remind the reader that expressions for the MBCS theory were copied from zero-temperature 
BCS expressions by straightforward replacing the BCS {uj,Vj} coefficients with the MBCS 
{uj, Vj} coefficients. It has been already pointed out [BJ that a method of mechanical copying of 
equations from one theory into another one has absolutely no grounds because the properties of 
{uj,Vj} and {uj,Vj} coefficients are simply different. E.g., the thermal part of the pairing gap 
5 A, discussed in Sec. El has appeared in the MBCS theory as a result of such an uncontrolled 
copying. 

It is very natural to verify whether expressions for thermodynamical observables obtained by 
such a method are thermodynamically consistent in the MBCS theory. Fig. [6] compares the 
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Fig. 6. MBCS predictions for the system entropy calculated as 5th, Eq. ((TJ), and S sp , Eq. ([8]) in different 
light and heavy, symmetric and asymmetric PFM systems. Notice the logarithmic y scale. 

system entropy which is calculated according its thermodynamical definition (dashed lines): 

J 1 dE tot 

Sth = J r'^r dT 

o 

and within a statistical approach (solid lines): 

S sp = -J2^J + l)[n j \nn, + (l-n 3 )ln(l-n j )] . (8) 

3 

The latter one is also called the single particle entropy in Refs. [TOlfTT] or the quasiparticle 
entropy in Ref. [8]. 

Fig. [6] shows that S t h and S sp being considered separately are amazingly close in different PFM 
systems for T < 1 MeV but they are dramatically different one from another in the MBCS 
predictions. The latter fact has been already addressed [5] when S"^ 803 and S^ 603 for the 
neutron system of 120 Sn have been found different by two orders of magnitude. 

Different entropy-like quantities in nuclear physics have been considered in [TOfTTj . It has been 
found that "the thermodynamic entropy . . ., the information entropy . . . and the single-particle 
entropy . . ., all coincide for strong enough interaction but only in the presence of a mean field" 
[TTj . It is not possible to accept interpretation of these results in Ref. [5] as that and S ap 
"are nearly the same only for noninteracting particles". Notice, correspondence between 
and S sp in a mean field plus residual forces (Figs. 56(IIa, lie) in [TU]) is definitely not worse 
than in an almost pure mean field approach (Figs. 56(Ia, Ic) in [TU]) . 

Fig. [6](a) also presents the quantities S^p s and Sf p cs calculated within the conventional BCS 
theory and S^ act obtained from the exact solution of the PFM for the N = Q = 10 system. 
One notices that for T > 100 keV: 

C-exact ^ qBCS ^ r-BCS ^ qMBCS r-MBCS 
°th ~ °th ~ °sp ~ °sp ^ ^ °th 




(7) 



(9) 



An attempt to compare S^ 803 and S^ 808 with the entropy calculated from the exact solution 
of the PFM has been made in Fig. 6 of Ref. [8] for the N = 10, Q = 11 system. First of all, the 
calculation of 5|p act (thick solid line in Fig. 6(a) of Ref. [8]) is definitely not correct. The author 
of [8] claims that S^ act does not vanish at T = because occupation probabilities fh < 1 and 
f p > 0. The problem with the third law of thermodynamiclf] in this calculation is caused by 
confusion of interacting particles with noninteracting "quasiparticles" : particle levels are not 
eigen states of the pairing Hamiltonian, their occupation numbers do not obey Fermi-Dirac 
distribution and because of that, they cannot be used in Eq. ([8]) which represents the free 
Fermi-gas combinatorics. Incorrectly calculated S^ act is also published in Fig. 8 of [3]. 

One notices that Fig. 6(b) of Ref. [S] lacks comparison between >S' t I ^ BCS and 5*™°* quantified" 3 !. 
The latter quantity can be easily calculated from E^ act (thick solid line in Fig. 4(a) of Ref. [8]). 
A level of expected disagreement between 5' t I ^ BCS and S^ act quantities can be seen from our 
Fig. El(a) (compare dashed curve with circles). 

Analysis of the MBCS predictions for the system entropy clearly indicate a problem with the 
expression for the system energy which enters in Eq. ([7]). How this expression was obtained, is 
mentioned in the beginning of the section. In Fig. 9 of Ref. [6] the system energy calculated 
from this expression has been straightforwardly compared to the system energy calculated as 
< H >= Tr(HD) where D is a density operator and < . . . > means averaging over the grand 
canonical ensemble. Dramatic disagreement in two quantities representing the same physical 
observable has been obtained indicating the same problem. 



6 Conclusions 



In this article we continue discussion on the validity of the MBCS theory. The theory perfor- 
mance is examined within the PFM. Our present goal is to allow the reader, who is not familiar 
with the previous discussion on the subject in Refs. [lfo]l6]r7] . to better understand the results 
and conclusions in Refs. [3"f8"] . 

We confirm that there exists a single example of the PFM (with the number of levels Q equal to 
the number of particles iV plus one) in which the MBCS produces the thermal behavior of the 
pairing gap similar to the one of a macroscopic theory up to a rather high temperature. In our 
opinion, this fact alone is not sufficient to conclude on the theory performance without noticing 
that in all other examples of the PFM with Q ^ N the theory predicts phase transitions of 
unknown types at a much lower temperature. These other examples are neither shown nor 
discussed in Ref. [3]. 

Although the MBCS formally yields a visually acceptable pairing gap behavior in some tem- 
perature domain, the physical content of its pairing gap is very dubious. The MBCS predicts 
that heating generates a thermal constituent of the pairing gap which becomes stronger and 
stronger with temperature in contradiction with generally accepted understanding of the pair- 
ing phenomenon in nuclei as a result of a specific particle-particle interaction. We point out 
that this unphysical constituent of the MBCS pairing gap is responsible for strong sensitivity 



"at absolute zero, any part of the body must be in a definite quantum state - namely the ground 
state . . . the entropy of the body - the logarithm of its statistical weight - is equal to zero" (see §23, 
page 66 in [12]). 

3 What is plotted by thick solid line in this figure and makes impression of agreement with S*^ BCS is 
not specified. 



of the theory predictions to tiny details of single particle spectra. 

The conventional PFM with Q = N possesses internal particle-hole symmetry. Due to it, some 
quantities in the model are known exactly without any calculations. We have demonstrated 
that the MBCS predictions deviate from exact results by a few hundred per cent starting from 
very low temperatures. 

We also point out that the MBCS is a thermodynamically inconsistent theory. 

Performing a systematic analysis of the MBCS performance in the PFM we have failed to find 
a model system and/or temperature range where the MBCS predictions are not questionable. 
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